import numpy as np
import matplotlib.pyplot as plt
# 定义函数f(k, q)
def f(k, q):
    return 1 + (1 / k) - q ** k
# 定义p值(阳性概率)
p = 0.0001
# 阴性概率
q = 1 - p
# 定义x的取值范围
k = np.arange(0,301)
y = f(k, q)
plt.plot(k, y)  # 添加数据
# 设置x轴的刻度间隔为5
plt.xticks(np.arange(0, 301, 50))
# 设置标题
plt.title('p=0.0001')
# 显示图形
plt.show()